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Abstract 

We calculate the orbital magnetization of a confined 2DEG as a function of the 
number of electrons in the system. Size effects are investigated by systematically 
increasing the area of the confining region. The results for the finite system are 
compared to an infinite one, where the magnetization is calculated in the thermo- 
dynamic limit. In all calculations the electron-electron interaction is included in the 
Hartree approximation. 
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1 Introduction 



Magnetization measurements offer a way to directly probe the equilibrium 
ground state properties of 2DEGs in the quantum Hall regime. Such informa- 
tion is not directly accessible in magnetotransport or optical spectroscopy, in 
which non-equilibrium behavior of the system is measured. Many years ago the 
de Haas- van Alphen oscillations of the magnetization were observed and used 
to characterize the DOS of the 2DEG [1]. In recent years many-body effects, 
including those related to the FQHE and skyrmions, have also been observed 
in experiments [2,3]. Two different experimental setups are used: one method 
uses sensitive mechanical, torque magnetometers [1,2], and a more recent one 
uses a low noise, superconducting quantum interference device (SQUID) [3,4]. 
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Various thermodynamic quantities, including the magnetization, have been 
calculated both for a disordered homogeneous 2DEG within the Hartree-Fock 
approximation (HFA) [5] and a statistical model for inhomogeneities corre- 
sponding to the Hartree approximation (HA) [6]. These calculations show 
many-body effects reflecting screening and, in addition, the HFA produces ex- 
change enhancement at odd filling factors due to spin splitting, and at even 
filling factors an enhancement, explored to lesser extent in experiments. 

For finite systems Kotlyar et al. calculated the persistent current in an array 
of quantum dots using Mott-Hubbard model for the intra-dot and inter-dot 
electron interaction [7,8]. For two spin split levels per dot oscillations in the 
persistent current, which is equivalent to magnetization density, were observed 
as a function of the number of flux quanta penetrating the system. 

In the present work we calculate the magnetization of a two dimensional 
electron system confined to a rectangular box, subjected to a finite one- 
dimensional periodic potential. Starting from a single unit cell we progressively 
increase the number of cells until we find convergence of the magnetization to 
the saw-tooth profile expected from the thermodynamic results. We discuss 
the edge and the bulk contribution to the orbital magnetization and the effects 
of the potential. We calculate separately the effects of the periodic potential in 
the infinite system, and compare the results. The electron-electron interaction 
is calculated in the HA for both the finite and the extended system. 



2 Models 

The model for a finite system consists of a laterally confined 2DEG. A hard 
wall potential ensures that the electrons stay in the rectangular region 



A perpendicular magnetic field and an external modulating potential are ap- 
plied to the system. The potential, which is of the form 



models Ux {rix integer) parallel quantum wires of width — Lx/rix and length 
Ly. The length is chosen to be Ly = UyC.j. so the geometry of system is con- 
trolled by the parameters and ny. By defining a unit cell of area the 
parameters Ux and Uy count the number of cells in the x and y direction, 
respectively. The total number of cells is N^y — UxUy. The wave functions of 



E = {(x,7/)|0 <x < Lx,Q <y < Ly}. 



(1) 
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the finite systems are expanded in sine Fourier series, since they are required 
to vanish at the boundary and the Schrodinger equation is solved in the HA. 
The orbital magnetization in the finite system is calculated according to the 
definition [9] 

M = -1- / d^r(r X (J(r))) • e,. (3) 



In the extended system the ground state is calculated in the HA, by diago- 
nalizing the Hamiltonian in the Landau basis, and by expanding the matrix 
elements as Fourier series. Therefore we can evaluate directly the magnetiza- 
tion by the thermodynamic formula appropriate for the canonical ensemble, 

M^--^^{E-TS\ (4) 



where E is the total energy, and S the entropy. We shall assume that the 
temperature is sufficiently low to neglect the entropy term. 



3 Results 



We use GaAs parameters, m* — 0.067mo, n — VIA. The magnetic field is 
B — 1.5 T, resulting in a magnetic length — 21 nm. The width of the 
wires is ^x = 75 nm, modulation amplitude Vq = lOmeV, and the temperature 
T = 1 K. In order to compare the results for the finite system to the infinite 
one we define the quantity v = Ns/Nlb, where Ng is the number of electrons 
and Nlb = LxLy/2Tril is the Landau level degeneracy. When -C L^Ly we 
can interprete P as the filling factor u. 

Since the modulation potential is anisotropic the system will behave differently 
depending on whether Ux or Uy is increased. Here we consider the cases of 
one, two , and three parallel wires, rix — 1, 2 and 3, of varying length. We 
calculate for lengths, Ly — Uyix, corresponding to % = 1, 2, 3, and 4. The 
magnetization for the finite system is shown in figure 1. It is plotted in units 
of NxyMo, where Mq = fJ^^/^l and fi*^ is the Bohr magneton containing the 
effective electron mass. The top panel shows the magnetization of a single 
wire, Ux — 1- For few electrons, P < 2, the current density circulates clockwise 
around the maximum density of the wire but for P between 2 and 3, depending 
on ny, the direction of the current reverses and the magnetization decreases 
again. Around z/ = 4 a single, dominant current loop splits into two smaller 
ones for Uy > 1 (the current simply reverses for Hy = 1) and the magnetization 
increases. The magnetization of two parallel wires, Ux — 2, is shown in the 
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Fig. 1. The magnetization for different x Uy. In the top panel the magnetization 
for a single wire, = 1, of length Ix through 4 x is plotted. The center and 
bottom panel show similarly the magnetization for Ux = 2 and rix = 3, respectively. 

center panel in figure 1. More pronounced jumps around u = 2 and 4 are 
due to denser energy levels because of the increased system size. In figure 2 
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Fig. 2. Energy spectrum as a function of Ns/N^b for two wires n^; = 2 of length 
Hy = 4. The solid line denotes the chemical potential fi. 
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the energy levels of the system and the chemical potential, shown as a solid 
line, are plotted as a function of the number of electrons. Below v ~ 2.5 the 
chemical potential traverses through relatively sparse (longitudinal) energy 
levels until it enters an energy interval corresponding to the second transverse 
energy band, with a lower slope. Around z/ = 4 it grows again. Screening 
reduces the amplitude of the modulation potential, but for < 3 the wires 
are well isolated. A more complex situation occurs for three wires, = 3, 
since there is a wire in the center of the system which is weaker coupled 
to the edges than the other two. At the boundary of the system there is a 
minimum in the Hartree potential due to the positive background charges 
[10]. As the wire length is increased, more states get localized around this 
minima and their contribution to the magnetization becomes more important. 
An interplay of this effect and contribution from the center wire causes the 
magnetization minimum preceeding the jump at = 4 to be shifted upwards, 
and for = 4 the region between u = 2 and 4 is relatively flat. In order 
to estimate the bulk contribution to the magnetization we use equation (3), 
but instead of integrating over the whole system we define a bulk area which 
is intergrated over. We choose this area to cover the center wire where we 
have taken 40 nm of its ends. The bulk magnetization, using the previous 
definition, is shown in figure 3. This can be compared with the magnetization 
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Fig. 3. The bulk magnetization for tree parallel wires of length Uy = 3 and Uy = 4. 
See text for the definition of the bulk contribution. 

calculated for an infinite system of parallel wires (n^ = oo, Uy = oo) also 
defined by the potential (2), shown in figure 4. In this case the energy spectrum 
consists of periodic Landau bands. The deviations from the saw-tooth profile, 
corresponding to the homogeneous system (Knod = 0) are determined both 
by the reduced energy gaps, i. e. the reduced jumps, and also by the energy 
dispersion, i. e. the nonlinear behavior in between even-integer z/-values. 

In summary, the wire modulation used here shows that magnetization due 
to edge states is nontrivial and large in mircroscopic systems. At the same 
time the bulk contribution to the magnetization for the system assumes very 
early the form known for an infinite 2DEG when the system size is increased. 
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Fig. 4. Thermodynamic magnetization for a modulation with Ix = 75 nm, in the 
HA. The solid line is Vq = lOmeV and dashed line Vq = 5meV. 

This might suggest experiments with a SQUID loop placed inside the system 
boundaries in order to verify the different contributions to the magnetization. 
The measurements of magnetization seem to provide a more direct access to 
the intrinsic equilibrium quantum many-electron structure of the system than 
transport and far-infrared experiments. 
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